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Abstract. Limiting cases are studied of the Koornwinder-Macdonald multivari- 
able generalization of the Askey- Wilson polynomials. We recover recently and not 
so recently introduced families of hypergeometric orthogonal polynomials in several 
variables consisting of multivariable Wilson, continuous Hahn and Jacobi type poly- 
nomials, respectively. For each class of polynomials we provide systems of difference 
CO , (or differential) equations, recurrence relations, and expressions for the norms of the 

polynomials in terms of the norm of the constant polynomial. 
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1. Introduction 

It is to date over a decade ago that Askey and Wilson released their famous memoir 
O ! [|AW|] , in which they introduced a four-parameter family of basic hypergeometric 

13 



polynomials nowadays commonly referred to as the Askey- Wilson polynomials ||GR | 



These polynomials, which are defined explicitly in terms of a terminating 403 series, 
have been shown to exhibit a number of interesting properties. Among other things, 
it was demonstrated that they satisfy a second order difference equation, a three- 
term recurrence relation, and that — in a suitable parameter regime — they constitute 
^ ■ an orthogonal system with respect to an explicitly given positive weight function with 

^ ! support on a finite interval (or on the unit circle, depending on how the coordinates 

^ I are chosen). 

Many (basic) hypergeometric orthogonal polynomials studied in the literature arise 
as special (limiting) cases of the Askey- Wilson polynomials and have been collected in 
the so-called (g-) Askey scheme [ |AW| , [KS|| . For instance, if the step size parameter of 
the difference equation is scaled to zero, then the Askey- Wilson polynomials go over 
in Jacobi polynomials: well-known classical hypergeometric orthogonal polynomials 
satisfying a second order differential equation instead of a difference equation. One 
may also consider the transition from orthogonal polynomials on a finite interval to 
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orthogonal polynomials on a (semi-)infinite interval. This way one arrives at Wilson 
polynomials (semi-infinite interval) and at continuous Hahn polynomials (infinite 
interval). 

The purpose of the present paper is to generalize this state of affairs from one to 
several variables. Starting point is a recently introduced multivariable generalization 
of the Askey- Wilson polynomials, found for special parameters by Macdonald | [[V12|| 
and in full generality (involving five parameters) by Koornwinder By means of 
limiting transitions similar to those in the one- variable case, we arrive at multivariable 
Jacobi polynomials |V|, pe|| (see also |[BU|| and reference therein) and at multivariable 
Wilson and continuous Hahn polynomials |p3| 



The (g-)Askey scheme involves many more limits and special cases of the Askey- 
Wilson polynomials than those described above. For instance, one also considers 
transitions from certain polynomials in the scheme to similar polynomials with less 
parameters and transitions from polynomials with a continuous orthogonality mea- 
sure to polynomials with a discrete orthogonality measure. Such transitions (or rather 



their multivariable analogues) will not be considered here. We refer instead to ||DT 
Sec. 5.2] for the transition from multivariable Askey- Wilson polynomials to Mac- 
donald's g-Jack polynomials (i.e., multivariable g-ultraspherical polynomials) |[M4|| 
(as an example of a limit leading to similar polynomials but with less parameters), 
and to ||SK|| for the transition from multivariable Askey- Wilson polynomials to mul- 
tivariable big and little g- Jacobi polynomials [|^] (as an example of a limit leading to 
multivariable polynomials with a discrete orthogonality measure). 

Whenever one is dealing with orthogonal polynomials an important question arises 
as to the explicit computation of the normalization constants converting the poly- 
nomials into an orthonormal system. For Jacobi polynomials calculating the or- 
thonormalization constants boils down to the evaluation of (standard) beta integrals, 
whereas Askey- Wilson polynomials give rise to g-beta integrals. In the case of sev- 
eral variables one has to deal with Selberg type integrals (Jacobi case) and g-Selberg 
type integrals (Askey- Wilson case), respectively. For these multiple integrals explicit 
evaluations have been conjectured by Macdonald that were recently checked using 
techniques involving so-called shift operators [Pp| , pS| , [Clj |M5|| . (Roughly speaking 
these shift operators allow one to relate the values of the (g-) Selberg integral for 
different values of the parameters separated by unit shifts; the integral can then be 
solved, first for nonnegative integer-valued parameters by shifting the parameters to 
zero in which case the integrand becomes trivial, and then for arbitrary nonnegative 
parameters using an analyticity argument (viz. Carlson's theorem).) 

Very recently, the author observed that Koornwinder's second order difference 
equation for the multivariable Askey- Wilson polynomials may be extended to a sys- 
tem of n (= number of variables) independent difference equations |pi|| and that the 
polynomials also satisfy a system of n independent recurrence relations ||D5||. (To date 
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a complete proof for these recurrence relations is only available for a self-dual four- 
parameter subfamily of the five-parameter multivariable Askey- Wilson polynomials.) 
It turns out that the recurrence relations, combined with the known evaluation for 



the norm of the unit polynomial (i.e., the constant term integral) [pu| , [Ka]| , may 
also be used to verify Macdonald's formulas for the orthonormalization constants of 
the multivariable Askey- Wilson polynomials | |D5| | . Below, we will use these results 
to arrive at systems of difference (or differential) equations, recurrence relations and 
expressions for the orthonormalization constants, for all three limiting cases of the 
multivariable Askey- Wilson polynomials considered in this paper (Wilson, continuous 
Hahn and Jacobi type). 

We would like to emphasize that much of the presented material admits a phys- 
ical interpretation in terms of Calogero-Sutherland type exactly solvable quantum 
n-particle models related to classical root systems ||0P|| or their Ruijsenaars type dif- 
ference versions |[R1| , [R2| , P2|| . The point is that the second order differential equation 
for the multivariable Jacobi polynomials may be seen as the eigenvalue equation for a 
trigonometric quantum Calogero-Sutherland system related to the root system BCn 
P| . From this viewpoint the second order difference equation for the multivariable 



Askey- Wilson polynomials corresponds to the eigenvalue equation for a Ruijsenaars 
type difference version of the i^C^-type quantum Calogero-Sutherland system [p4 |. 
The transitions to the multivariable continuous Hahn and Wilson polynomials amount 
to rational limits leading to (the eigenfunctions of) similar difference versions of the 
A„_i-type rational Calogero model with harmonic term (continuous Hahn case) and 
its B{C)n-tYpe counterpart (Wilson case) |p3|| . For further details regarding these 
connections with the Calogero-Sutherland and Ruijsenaars type quantum integrable 
n-particle systems the reader is referred to [P^ , P^ , |D4|| . 

The material is organized as follows. First we define our families of multivari- 
able (basic) hypergeometric polynomials in Section ^ and recall their second order 
difference equation (Askey- Wilson, Wilson, continuous Hahn type) or second order 
differential equation (Jacobi type) in Section |^. Next, in Section ^, limit transitions 
from the Askey- Wilson type family to the Wilson, continuous Hahn and Jacobi type 
families are discussed. We study the behavior of our recently introduced systems of 
difference equations and recurrence relations for the multivariable Askey- Wilson type 
polynomials with respect to these limits in Sections ^ and ^ respectively. The recur- 
rence relations for the Wilson, continuous Hahn and Jacobi type polynomials thus 
obtained in Section ^ are then employed in Section ^ to derive explicit expressions 
for the (squared) norms of the corresponding polynomials in terms of the (squared) 
norm of the unit polynomial. 
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2. MULTIVARIABLE (bASIC) HYPERGEOMETRIC POLYNOMIALS 

In this section multivariable versions of some ortliogonal families of (basic) hyper- 
geometric polynomials are characterized. The general idea of the construction (which 
is standard, see e.g. |V|, |M2| , |K|, \SK\\ ) is to start with an algebra of (symmetric) poly- 
nomials Ti. spanned by a basis of (symmetric) monomials {'mx}x^\, with the set A 
labeling the basis elements being partially ordered in such a way that for all A G A the 
subspaces Hx = Span{m^}^gA,/i<A are finite-dimensional. It is furthermore assumed 
that the space TC is endowed with an inner product (-, ■)a characterized by a cer- 
tain weight function A. To such a configuration we associate a basis {px}x£A of Ti 
consisting of the polynomials px-, A G A, determined (uniquely) by the two conditions 

I. px = mx+ ^ Ca,^ m^, cx,^, G C; 

fj.i=A,ii<X 

a. {px,mij)A = if /i < A. 

In other words, the polynomial px consists of the monomial rrix minus its orthogonal 
projection with respect to the inner product (■, ■)a onto the finite-dimensional sub- 
space Span{m^}^gA,/x<A- By varying the concrete choices for the space H, the basis 
{mA}AeA and the inner product (■,-)a, we recover certain (previously introduced) 
multivariable generalizations of the Askey- Wilson, Wilson, continuous Hahn and Ja- 
cobi polynomials, respectively. Below we will specify the relevant data determining 
these families. The fact that in the case of one variable the corresponding polynomi- 
als Pa indeed reduce to the well-known one-variable polynomials studied extensively 
in the literature is immediate from the weight function. The normalization for the 
polynomials is determined by the fact that (by definition) px is monic in the sense 
that the coefficient of the leading monomial nix in Px is equal to one. 

It turns out that in all of our cases the basis {mAjAgA can be conveniently expressed 
in terms of the monomial symmetric functions 

msymAzi, . . . , Zn) = Yl ^ ^ ^' (2-1) 

fieSn{X) 

where 

A = {A G Z" I Ai > A2 > • • ■ > A„ > 0}. (2.2) 
In (|2.1|) the summation is meant over the orbit of A under the action of the permuta- 



tion group Sn (acting on the vector components Ai, . . . , A^). As partial order of the 
integral cone A ( |2.2| ) we will always take the dominance order defined by 

m m 

fi < X iff fij < Xj for m = 1, . . . ,n (2-3) 

i=i i=i 

(and /i < A iff /i < A and /i 7^ A). 
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Note: In order to avoid confusion between the various families we will often equip 
the polynomials and other objects of interest with the superscripts 'AW, 'W, 'cH' or 
'J' to indicate Askey- Wilson, Wilson, continuous Hahn or Jacobi type polynomials, 
respectively. Sometimes, however, these superscripts will be suppressed when dis- 
cussing more general properties of the polynomials that hold simultaneously for all 
families. 



2.1. Askey-Wilson type. To arrive at multivariable Askey- Wilson type polynomi- 
als one considers a space TC"^^ consisting of even and permutation invariant trigono- 
metric polynomials. Specifically, the space TC^^ is spanned by the monomials 

mf^{x) = m,j,„,A(e'""^ + e-'""\ • • • , e'""" + e"'"""), A G A (2.4) 

(with A given by ( |2.2| )). The relevant inner product on T-C^^ is determined by 

J-TT/a J-n/a (2.5) 



with the weight function reading 



; . , ^ ' (2.6) 



l<j<k<n 
ei,£2 = ±l 

^ i-r (e^-"^"^;g)oo 

J- J- (to e'°^^J , ti e^'^'^^i , t2 e^'^^'^i , ts e*"=^i ; q) c 

l<j<n 

e=±l 



Here (a; g)oo = 11^=0(1 "^'")' ("i' • • • ^^r] q)oo = (en; g)oo ■ ■ ■ (ctr; g)oo and the pa- 
rameters are assumed to satisfy the constraints 

a > 0, 0<g<l, -l<t<l, < 1 (r = 0, 1,2,3), (2.7) 

with possible non-real parameters tr occurring in complex conjugate pairs and pair- 
wise products of the tr being ^ 1. For the weight function in (|2.6|) the polynomials 
Px determined by the Conditions i. and ii. (above) were introduced by Macdonald 
[|M2|| (for special parameters) and Koornwinder [Q (for general parameters). In the 
special case of one variable {n = 1) these polynomials reduce to monic Askey-Wilson 
polynomials ||AW| , |KS|| 



AWf \ _ (^o^i! ^0^2, ^0^3! Q')i , ((1 \ tQtit2t-iq^ ^, to^*"^, t^e ^ 

Pi ~ Trmm. — m — r 4</^3 ;9'9' 



to{totit2t3q^-^;q)i \ toti, tot2, tot 
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2.2. Wilson type. In the Wilson case the appropriate space Ti.^ consists of even 
and permutation invariant polynomials and is spanned by the monomials 

mY{x) = nisyraAxl,- ■ ■ ,xl), A G A. (2.9) 
The inner product on Ti^ is now determined by 



(mf , m^)^w =/■■■/ mY{x) mjf (x) A^(x) dxi--- dxn, (2.10) 

with the weight function taken to be of the form 
A^(x) = n n- + ^{e^x,+e,x,)) 



l<j<k<n ^ ^ * J 



T{i{eiXj + e2Xk)) 

ei,e2=±l 

r(z/o + ieXj)T{h'i + ieXj)T{h'2 + iexj)T{h'3 + iexj) 



n 



T(2ieXi 

l<j<n ^ •> 

e=±l 



Here r(-) denotes the gamma function and the parameters are such that 

z/>0, Re(i/^) > (r = 0,1,2,3), (2.12) 

with possible non-real parameters Vr occurring in complex conjugate pairs. For this 
weight function (and real parameters) the polynomials p\ were introduced in [D3|. 
In the case of one variable they reduce to monic Wilson polynomials [|AW| , |KS 



4-^3 



z/q + z/i + z/2 + z/3 + / - 1, z/q + ix, VQ-ix _^ 

^0 + Z^l, + ^2, Vq + V^ 



2.3. Continuous Hahn type. The space T-f^^ is very similar to that of the Wilson 
case but instead of only the even sector it now consists of all permutation invariant 
polynomials. The monomial basis for the space l-f^^ then becomes 

mf^{x) =msymA^ir-- ,Xn), A G A. (2.14) 
The inner product is of the same form as for the Wilson case 



(m^^, mf)^.H =/■■•/ mf{x) mf{x) A^^(a;) rfxi ■ ■ ■ dxr. 

but now with a weight function given by 

.cH/^^ _ jj- fT{u + i{xj - Xk))T{iy + i{xk - Xj)) 

l<j<k<n 

i<i<" 



(2.15) 



' 11 I Y{i{x^-Xk)) V{i{xk-x,)) I ^ ■ ^ 
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where 



u>0, Re(z/o^),Re(z/i^) > 0, 



'0 ' 



(2.17) 



Just as in the case of Wilson type polynomials, the polynomials corresponding to 
the weight function ( p.l6|) were introduced in |P3|| . For n = 1 they reduce to monic 
continuous Hahn polynomials |[A"W| , [KS|| 



i' ji^^ + Up , t^o+ + )i 



X 



(2.18) 



3-^2 



' 



Z/i 



2.4. Jacobi type. The space Ti^ and the basis {fri'l^Xf^K are the same as for the 
Askey- Wilson type. The inner product is also of the form given there (cf. ( p.5|) ) but 
the weight function gets replaced by 



n 

l<j<k<n 

- n 

i<i<™ 



sm-(Xj +a;fc)sm-(Xj 



sin(-a;j, 



21^0 



cos(-Xj, 



2vi 



2v 



(2.19) 



with 



a > 0, z/>0, z/o,z^i > -1/2. (2.20) 

In this case the corresponding polynomials p\ were first introduced by Vretare |V| 
(see also |P4 [B(J|| ). For n = 1 they reduce to monic Jacobi polynomials ||A!j| , [KS|| 

2^' (//o + 1/2), 



2^1 



-/, Z/Q + Z/i + / .2 

z/o + 1/2 ' '''' V 2 y 



(2.21) 



(z/Q + Z/1 + /); 

Remark: In the explicit formulas for the polynomials when n = 1 we have used 
their standard representations in terms of terminating (basic) hypergeometric series 



AW, KS 



F 

r-'- s 



s+l<+>s 



hi,... ,b, 

Oi, . . . , Og+i 

bi,...,bs ' 



; z 



q,z 



E 

k=0 

oo 

E 

A:=0 



(ai, . . . ,ar)k Z^ 
{bi, . . . ,bs)k k\ ' 

(qi, . . . , g^+i; q)k z^ 
(6i,... ,bs]q)k {q;q)i 



(2.22) 
(2.23) 



where we have used Pochhammer symbols and g-shifted factorials defined by 

(ai, . . . , ar)k = (ai)fc ■ ■ ■ (ar)fc, (^i, • • • ,ar] q)k = [ai] q)k-- ■ [ar] q)k 
with (a)o = (a; q)o = 1 and 

{a)k = a{a + 1) ■ ■ ■ {a + k - 1), (a; q)k = (1 - a)(l - ag) • ■ ■ (1 - ag''""^) 
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for A; = 1, 2, 3, ... . 



3. Second order difference or differential equations 

As it turns out, all families of polynomials {px}x<=a introduced in the previous 
section satisfy an eigenvalue equation of the form 

Dpx = ExPx, AeA, (3.1) 

where D -.Ti ^Ti denotes a certain second order difference operator (Askey- Wilson, 
Wilson and continuous Hahn case) or a second order differential operator (Jacobi 
case). Below we will list for each family the relevant operator D together with 
its eigenvalues i?A, A G A. In each case the proof that the polynomials px indeed 
satisfy the corresponding eigenvalue equations boils down to demonstrating that the 
operator D : TC ^ H maps the finite-dimensional subspaces Hx = Span{m^}^gA,^<A 
into themselves (triangularity) and that it is symmetric with respect to the inner 
product (■, ■)a- In other words, one has to show that 
Triangularity 

Dmx= with [DK^eC (3.2) 

/i€A,/i<A 



and that 
Symmetry 



{Dmx,m^)A = {mx,Dm^)A- (3.3) 



It is immediate from these two properties and the definition of the polynomial px that 
Dpx lies in Hx and is orthogonal with respect to (■, ■)a to all monomials m^, /i G A 
with fi < X. But then comparison with the definition of px shows that Dpx must be 
proportional to px, i.e., px is an eigenfunction of D. The corresponding eigenvalue 
Ex is obtained via an explicit computation of the diagonal matrix element [-D]a,a in 
Expansion (|3.2|) . 

For the Jacobi proof of the second order differential equation along the 

above lines was given by Vretare [0]. In the Askey- Wilson case the proof was given 



by Macdonald [M2] and (in general) Koornwinder [|^. The proof for the Wilson and 



continuous Hahn case is very similar to that of the Askey- Wilson case and has been 



outlined in D3 . 



3.1. Askey- Wilson type. The second order (g-) difference operator diagonalized by 
the polynomials p^^, A G A, is given by 

^Aw ^ ^ (l^^^(x)(T,,, - 1) + V^^{x){Tr^' - 1)) (3.4) 
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with 

^AW I \ _ ii0<r<3l 



n/1 _ J. ±iax 



X 

l<fc<n, fc^j 



Here the operators T,- ^ act on trigonometric polynomials by means of a g-shift of the 
jth variable 

(Tj J)(e^°^S . . . , e^°^") = /(e*"^'S . . . , e^"^'^-i, ge^"^% e^"^^+S . . . , e^°^"). 

(3.6) 

The eigenvalue of D"^^ on has the value 

Ex"^ = E {totit2hq~h'-~'~\q'^ - 1) + - 1)) . (3.7) 

i<j<" 

Proposition 3.1 ([Q). T/ie multivariahle Askey-Wilson polynomials Px^, A G A 
(p.2|) , satisfy the second order difference equation 

DAWpAw ^ E^'^pf^. (3.8) 



For n = 1, Equation (|3.8| ) reduces to the second order difference equation for the 
one-variable Askey-Wilson polynomials [|AW|, [KS| 



^ ^ (pr(gx)-pr(-)) + (3.9) 



1 - e2*°^) (1 - ge2»"^y 



-tax \ 



(1 - e-2»°^) (1 - ge-2«°^'; 

= (toW3g''(g' - 1) + (g-' - 1)) pr(x). 

3.2. Wilson type. In the case of Wilson type polynomials the difference operator 
takes the form 

D^=Y, [Vr^^m - 1) + V^^{x){T-' - 1)) (3.10) 



where 



X n 

l<fc<n, fc^i^j 



it/ ± -|- Xjt \ f iu ± Xj — Xk 



I 3/ ^ ^ \^ _l — 
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and the action of Tj is given by a unit shift of the jth variable along the imaginary 



axis 



(Tjf){xi, ... ,Xn) = /(Xi, . . . ,Xj-i,Xj + i,Xj+i, . . . ,Xn). (3.12) 

The corresponding eigenvalues now read 

= Yl i^J + Uo + i^i + i^2 + i^3~l + 2{n- j» . (3.13) 

l<i<n 

Proposition 3.2 ( [p3|| ). The multivariahle Wilson polynomials pY , A G A ( |2.2| ), 

satisfy the second order difference equation 

D^pY = EYpY. (3.14) 

For n = 1, Equation (|3.14|) reduces to the second order difference equation for the 
one-variable Wilson polynomials [[KS|] 

2ix(2.x-l) (^ + ^)-^^ (^)) + 



(3.15) 



(pj^(x-i) -^^(x)) 



2%x {2ix + 1) 
= / (Z + z/0 + z^i + z/2 + z/3 - 1)pY{x). 

3.3. Continuous Hahn type. For the continuous Hahn type one has 

D'"" = E i^l-i^m - 1) + ySi^)iTf' - 1)) (3.16) 
i<i<" 

with 

V^^ix) = {u^+zx,){ut + zx,) n fl + ,Y. ' (3.17) 

i<fc<",Mi ^ ^ ^' ''^^ 

\A^^(x) = (z.--za;,)(z/r-za;,) J] (^' u 'I " ^^.IS) 

The action of Tj is the same as in the Wilson case (cf. ( |3.12| )) and the eigenvalues 
are given by 

= E + "0 + + '^0 +1^1-1 + 2{n - j» . (3.19) 

l<i<n 

Proposition 3.3 ( [P3|| ). The multivariahle continuous Hahn polynomials p'^^ , A G A 
(p.2|) , satisfy the second order difference equation 

JJcHpCH ^ Efpf. (3.20) 
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For n = 1, Equation (|3.20| ) reduces to the second order difference equation for tlie 



one-variable continuous Halin polynomials ||KS|| 

(z/o" - ix){v^ - ix) {pf{x + - p^ix)) + (3.21) 

(z/0+ + tx)iu+ + IX) {pf{x - 2) - pf{x)) 



I {I + 4 + vt + ^^0 + ^1 - 1) Pi 



cH , 



X) 



3.4. Jacobi type. In the case of multivariable Jacobi type polynomials the operator 
D diagonalized by p\ is given by a second order differential operator of the form 



Y^d] -aY^(u, cot(^) - tan(^)) d, (3.22) 
(cot(-(x, + x,)){d, + d,) + cot(-(x, - x,)M - d,)) 



D' = - 



—au 

i<j<k<n 

where dj = d/dxj. The eigenvalue of D"^ on pi takes the value 



= 5Z A, (A^. + z/o + ^^1 + 2(r2 - j» . (3.23) 

l<i<n 



Proposition 3.4 ([|V|, ^e|). T/ie multivariable Jacobi polynomials p^, A G A (|2.2|), 

satisfy the second order differential equation 



D'pi = Eipi (3.24) 

For n = 1, Equation (|3.24|) reduces to the second order differential equation for 
the one- variable Jacobi polynomials [|AS| , |KS| 



^^(x) - « z/ocot(— )-z/itan(— ) -P(x) (3.25) 



dx"^^ ' \ ^ 2 ' ^ 2 'J dx 

= l{l + VQ + vi)pi{x). 

4. Limit transitions 

The operator D of the previous section can be used to arrive at the following useful 
representation for the polynomials p\ (cf. ||M2| , |D5| , |SK|| ) 



n #^) -- (") 



Indeed, it is not difficult to infer that the r.h.s. of (|4.1|) determines a polynomial 
satisfying the defining properties i. and ii. stated in Section |[ To this end one uses 
the Triangularity (|3.2|) and Symmetry ( p. 3D of D, together with the observation that 
in each of the concrete cases discussed above the denominators in (14.11) are nonzero 
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since for parameter values indicated in Section ||one has that (see |pi| , Sec. 5.2] and 

m) 

fi<X =^ E^<Ex. (4.2) 

(It is immediate from the triangularity of D that the r.h.s. of (^4.1|) can be written as a 
linear combination of monomials with /i < A; that the r.h.s. is also orthogonal to 
all rrin with /i < A follows from the symmetry of D and the fact that the operator in 
the numerator — viz. Yl^eA fj.<xi^ ~ -^m) — annihilates the subspace Span{m^}^gA^^<A 
in view of the Cayley-Hamilton theorem.) 



Below we will use Formula ( |4.1| ) to derive limit transitions from the Askey- Wilson 
tj^e to the Wilson, continuous Hahn and Jacobi type families, respectively. The 
transition 'Askey- Wilson Jacobi' has already been considered before in ||1V12|, [UT 



SK[| and is included here mainly for the sake of completeness. It will be put to use in 
Section ^ when deriving a system of recurrence relations for the multivariable Jacobi 
type polynomials. 

4.1. Askey- Wilson— > Wilson. When studying the limit — > it is convenient 
to first express the multivariable Askey- Wilson polynomials in terms of a slightly 
modified monomial basis consisting of the functions 



mf^(x) = (2/«)2|^l msy^,xism\axi/2), ... , sin2(ax„/2)), A G A, 



(4.3) 



where |A| = Ai + ■ ■ ■ + A„. Notice that 

\imm^^{x) =mY{x) (4.4) 

whereas the original monomials mf^^{x) ( p.4|) all reduce to constant functions in this 
limit. Using the relation sin^(Q!Xj/2) = 1/2 — (e*"'^-'' + e~*°^^)/4 one easily infers that 
the bases {mf^^}xi=A and {mf^^}xeA are related by a triangular transformation of 
the form 

rhx = (-l/a')l^l mf^ + ^ ax,^ with ax,^ G M. (4.5) 

It is clear that in Formula ( |4.1| ) we may always replace the monomial basis {^aIasa by 
a different basis that is related by a unitriangular transformation, since (cf. above) the 
operator YliieA fj.<x(^ ^ -^m) numerator of the r.h.s. annihilates the subspace 

Span{m^}^gA,fj<A because of the Cayley-Hamilton theorem. Hence, by taking in 
account the diagonal matrix elements in the basis transformation (|4.5| ), one sees that 
Formula ([4.1|) can rewritten in terms of as 

(fjAW _ pAW\ 
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If we now substitute 

g = e-", ^ = e-"^ = e""'^'' (r = 0, 1, 2, 3) (4.7) 
in D^^ (U) and then we have that 

hma-2D^^ = D^, lima-^E^"^ = EY. (4.8) 

a^O a— >0 

(Notice to this end that for q = e~" the action of Tj g (|3.6|) on trigonometric poly- 



nomials is the same as that of Tj ( |3.12|) , i.e., the action amounts to a shift of the 
variable xj over an imaginary unit: xj —>■ xj + i. To infer then that in the limit a —>■ 
the difference operator a~'^D^^ formally goes to boils down to checking that 
the coefficients of the operator converge as advertised.) 

By applying the limits ( [4. 41) and ( [4.8|) to Formula ( [4.6|) we end up with the fol- 
lowing limiting relation between the multivariable Askwey- Wilson and Wilson type 
polynomials. 

Proposition 4.1. For Askey-Wilson parameters given by ( [4.7| ) one has 

pY{x) = hm (-l/«2)l^l pr(^), A G A (4.9) 
(with |A| = Ai + ■ ■ ■ -|- Xn)- 

4.2. Askey-Wilson^continuous Hahn. Just like in the previous subsection, the 
derivation of the transition p'^^ hinges again on Formula ( [4.11) . If we shift the 

variables xi, . . . , x„ over a half period by setting 

Xj Xj — vr/ (2a), j = 1, . . . ,n (4-10) 

and substitute parameters in the following way 

q = e"", t = e-°^ 

U = -ie-'"'o, t^=-ie-'"'i, t-2 = ie~'^''o , = ie-°^i', 



(4.11) 



then the version of Formula ( |4.1| ) for the multivariable Askey-Wilson polynomials 
takes the form [cj denotes the jth unit vector in the standard basis of MJ^) 

^ ^ ^ga,m<a\^^ / (4.12) 

where 

^Aw ^ J2 (vf'^{x){T^-l) + r4^(x)(T-i-l)), 

i<i<" 
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with 



X n 

l<k<n,k=ij 



(1 + e2'"^i) (1 + e-"e2»"^0 



X n 

l<fc<n, fc^i^j 



X _|_ g-ia(xj+a;ft) / \ — g-«a(2;j -^^fc) 



and 



m^^{x) = nisymA'^ sin(axi), . . . , 2 sin(ax„)). 



(Just as in the case of the transition Askey- Wilson— Wilson we have rewritten the 
operators Tj^g (|3.6| ) for q = e~" as Tj (|3l2|) .) After dividing by (2a) 1^1 the r.h.s. 
of (|4.12|) goes for a — > to the corresponding formula for the continuous Hahn 
polynomials (i.e. with D^^ D"^ , Ef^ Ef and (2a)-l^lmf^ mf). 
Hence, we now arrive at the following limiting relation between the multivariable 
Askey- Wilson and continuous Hahn type polynomials. 

Proposition 4.2. For Askey-Wilson parameters given by ( [4.11|) one has 

where u; = ei + ■ ■ ■ + e„ (with Cj denoting the jth unit vector in the standard basis of 

4.3. Askey-Wilson^Jacobi. To recover the Jacobi type polynomials we substitute 
the Askey-Wilson parameters 

t = q^, to = q'\ t, = -q^\ t^ = q''^^'^\ t, = -q^'-^^/\ 

(4.14) 

With these parameters the formula of the Form ( |4.1| ) for the Askey-Wilson type poly- 
nomials reduces in the limit g — 1 to the corresponding formula for the Jacobi type 
polynomials (i.e., the difference operator D^"^ with eigenvalues E^"^ gets replaced 
by the differential operator D'^ with eigenvalues Ej^). The limit g — > 1 amounts 
to sending the difference step size to zero. In order to analyze the behavior of the 
operator D^"^ for g — 1 in detail it is convenient to substitute q = e~"^ (so the 
action of Tj g ( p.6|) on trigonometric polynomials amounts to the shift Xj Xj + i[3) 
and then write formally Tj^q = exp{i(3dj). A formal expansion in /3 then shows that 
jjAw ^ i^2jjj ^YisX Ef^ ~ 13'^Ei for /3 ^ 0. Here the Jacobi parameters i/, Vr 
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are related to the parameters g, gi'^ in ([4.14|) via = g, i^o = go + g'o and ui = gi + g'l 



As a consequence we obtain the following limiting relation between the multivariable 
Askey- Wilson and Jacobi type polynomials. 



Proposition 4.3. For Askey-Wilson parameters given by ( [4.14|) one has 

pi{x) = lim pf^(x) (4.15) 

with the Jacobi parameters v, Uq and Vi taking the value g, g^ + g^ and gi + g'l, 
respectively. 



Remarks: i. In the above derivations of Propositions ^]T|, and |4.3| we have 
used that the Askey-Wilson type difference operator converges formally (i.e., without 
specifying the domains of the operators of interest) to the corresponding operators 
connected with the Wilson, continuous Hahn and Jacobi type polynomials, respec- 
tively. In our case such formal limits get their precise meaning when being applied 
to Formula (|0|). 

a. For all our four families AW, W, cH and J the dependence of (the coefficients of) 
the operator D and of the eigenvalues Ex is polynomial in the parameters t, to, ti, t2, t^ 
(AW), z/, z/q, z/i, z/2, z/3 (W), z/, z/^,z/^ (cH) and z/, z/q, z^i (J), respectively. Hence it is 
clear from Formula ( [4.1[ ) that (the coefficients of) the polynomials p\ are rational in 
these parameters. We may thus extend the parameter domains for the polynomials 
given in Section ^to generic (complex) values by alternatively characterizing px as the 
polynomial of the form px = rnx + Xl^eA ijl<\ ^a.^j^^^i satisfying the eigenvalue equation 
Dpx = ExPx- It is clear that the limit transitions discussed in this section then extend 
to these larger parameter domains of generic (complex) parameter values. 

Hi. In the case of one variable the limit transitions from Askey-Wilson polynomials 
to Wilson, continuous Hahn and Jacobi polynomials were collected in | |KS[ (together 
with many other limits between the various (basic) hypergeometric orthogonal fami- 
lies appearing in the (g-) Askey scheme). 



5. Higher order difference or differential equations 

In [ pi[ | it was shown that the second order difference equation for the multivariable 
Askey-Wilson type polynomials can be extended to a system of difference equations 
having the structure of eigenvalue equations of the form 

DrPx = Er,xPx, r = l,...,n, (5.1) 

for n independent commuting difference operators -Di, -D2, • • • , of order 2,4,..., 2n, 
respectively. For r = 1 one recovers the second order difference equation discussed in 
Section |3.1| . After recalling the explicit expressions for the Askey-Wilson type differ- 
ence operators and their eigenvalues Ef^^ , we will apply the limit transitions 
of Section ^ to arrive at similar systems of difference equations for the multivariable 
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Wilson and continuous Hahn type polynomials. In case of the transition 'Askey- 
Wilson^ Jacobi' the step size is sent to zero and the system of difference equations 
degenerates to a system of hypergeometric differential equations, thus generalizing 
the state of affairs for the second order operator in the previous section. This limit 
from Askey- Wilson type difference equations to Jacobi type differential equations has 
already been discussed in detail in ||D1| , Sec. 4], so here we will merely state the results 
and refrain from presenting a complete treatment of this case. 

5.1. Askey- Wilson type. The difference operators diagonalized by the multivari- 
able Askey- Wilson polynomials via Eq. ( |5.1| ) are given by 

E UZ-\j\{^)yeZ^(^)TeJ,, r = l,...,n, (5.2) 

JC{1,... ,n}.0<\J\<r 

ej=±l, jeJ 

with 

^Jjjci^) = ^'"^'^'^(^i^i) Y\. '^'^^i^j^j + £j'Xj')v^'^ {SjXj +ej'Xj' — i\n{q)/a) 

j<j' 

je.J 

(-1)' E (n ^^"^(^i^i) n ^"""^i^i^i + ei'x,)v^'^{-eixi - evxv + i Hq)/a) 

LCK, \L\^p leL 1,1'eL 

X JJ v^"^ {eixi + Xk)v^'^ {eixi - Xk)^ , 

leL 
keK\L 

and 

Here the action of the operators T^^ is defined in accordance with ( |3.6| ). The sum- 
mation in (|5.2| ) is over all index sets J C {1, . . . , n} with cardinality | J| < r and over 
all configurations of signs Sj G {+1, —1} with j G J . Furthermore, by convention 
empty products are taken to be equal to one and Uk,p = 1 for p = 0. 
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The corresponding eigenvalue of on ^ has the value 



where 



and 



(5.5) 



Er{Xi, ... ,Xn;yr,--- .Vn) = (5.6) 

JC{1 n} ieJ f'<'l<---<'r-| /l^'^ 

0<|J1<T- 

= e~^{t^tihhq-Y\ J = 1, . . . , n. (5.7) 



(The second sum in ( |5.6|) is understood to be equal to 1 when \J\ = r.) 
Summarizing, we have the following theorem generalizing Proposition |3.1| . 



Theorem 5.1 ( |pi| , P4|| ). The multivariahle Askey-Wilson polynomials , A G A 
(p.2|) satisfy a system of difference equations of the form 



Drpr = E:iTpr, r = i,...,n. (5. 



For r = 1 the Difference equation ( p. 81 ) goes over in the second order Difference 
equation (|3.8| ) after multiplication by a constant with value t"'~^(to^i^2^3Q'~^)^^^- More 
generally, one may multiply the Difference equation ( |5.8| ) for arbitrary r by the con- 
stant factor t^("~^)~^*^''~^)/^(toiit2^3'?~^)^''^ to obtain a difference equation that is poly- 
nomial in the parameters t, t^, . . . , t^ and rational in q. Such multiplication amounts 
to omitting the factors and (to^i^2^3Q'^^)^^^^ in the definition of v^^{z) (|5.3|) 
and w'^'^(x) (|5^) and to replacing the eigenvalues by 



j^AW _^ ^-r(r-l)/2 



with 



xEr{r+q^' +T^q . . . , r+g^" + r'q r+ + r", . . . , r+ + r, 



r- = r-^(totit2t3g-')'/Vri = t^-i. 



5.2. Wilson type. If we substitute Askey-Wilson parameters of the form ( |4.7| ) and 
divide by a"^^, then for a — the operator D:^^ goes over in 

D^= E r = l,...,n, (5.9) 

JC{1,... ,n},0<|J|<r 
ej=±l, JG J 
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with 



j<j' 

X Y\ {SjXj + Xk)v^ {SjXj — Xk), 



jeJ 



[-ly ^ (^J]^ Y\ {^i^i + evXi>)v^ {-eiXi - evxv - i) 



LCK,\L\=p 1<^L l.l'eL 

ei=±l, l€L Kl' 



X 



Y[ v^{eixi + Xk)v^{eixi - Xk)h 



and 



keK\L 



Here the action of T is taken to be in accordance with (p.l2| ). To verify this hmit 



it suffices to recall that for g = e " the action of ^ ( |3.6| ) amounts to that of Tj 
(p.l2|) and to observe that for parameters ( [4.7|) one has \\ma->v^^{z) = v^{z) and 



lima^ a ^w^'^(z) = w^{z). 

The eigenvalues become in this limit 

<A = ((pf + x^)^ . . . , (pr + A.)^; (P^)^ ■ ■ ■ , (pD^) (s.h) 

with Er{- ■ ■ '■,■■■) taken from ( ^.6|) and 

Pf = - J> + ('^o + z/i + z/2 + ^^3 - l)/2. (5.12) 

For the eigenvalues the computation verifying the limit is a bit more complicated 
than for the difference operators; it hinges on the following lemma. 

Lemma 5.2 ([0, Sec. 4.2]). One has 



2 . „.2 



J2jr{Xi, . . . , X^, y^, . . . ^Un)- 

We may thus conclude that the transition AW^W gives rise to the following 



generalization of Proposition 3.2 
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Theorem 5.3. The multivariahle Wilson polynomials , A G A (|2^ ) satisfy a sys- 
tem of difference equations of the form 

D^pY = E^,pY, r = l,...,n. (5.13) 

For r = 1 the Difference equation ( |5.13| ) coincides with the second order Difference 
equation (p.l4|) . 



5.3. Continuous Hahn type. After shifting over a half period as in ([4.10|) and 
choosing Askey- Wilson parameters of the form ([4.11|) , the operators a~'^^D^^ go for 
a ^ over in 

= ^ UjfnJt,r-\j+\-\j^\i^)'^j+,j.-,Jinj±i^)'^J+,J- 

J+,J_C{1,... ,n} (5-14) 

J+nJ-=0, \ J+\ + \J-\<r 

r = 1, . . . ,n, where 

= n ^ ' n 



X 



kGK k€K 



Uf,,{x) 



L^,L_CK, L^nL_=li 1&L+ IgL- ;eL_|_ 

\L+\ + \L-\=p l'€L^ 

X n ^"""(^I'^k) n v'^^'i^k-xi)'), 



ieL_|_ ieL_ 
k&K\L+UL_ keK\L+UL. 



and 



v''"{z) = (1 + -), (5.15) 

wf{z) = {u+ + iz){u+ + iz), w'L^iz) = (z/q- - iz)(z/f - iz). (5.16) 

(The action of the operators is again in accordance with (|3.12|) .) To verify this 
transition one uses that v'^^{±{xj + Xk)) —>■ 1, ^^^{xj — Xk) — * v^^{xk — Xj) and 
that a~'^w'^^ {±Xj) — >• w^{xj), if one sends a to zero after having substituted ( ^4.10| ) 
and ( pll) . 

The computation of the limit of the eigenvalues is exactly the same as in the Wilson 
case and the result reads 

^3 = Er {{pf + Al)^ . . . , (pf + Xnf; {pf)\ (pf )') (5.17) 
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with Er{- ■ ■ ',■■■) taken from ( ^.61 ) and 

pf = in- J> + K + + '^0 + ^1 - l)/2- (5.18) 
Hence, we arrive the following generalization of Proposition |3.3| . 

Theorem 5.4. The multivariahle continuous Hahn polynomials , A G A (|2.2| ) 

satisfy a system of difference equations of the form 

Dfpf = Efpf, r = l,...,n. (5.19) 

For r = 1 the Difference equation (|5.19| ) coincides with the second order Difference 
equation (|3.20|) . 



5.4. Jacobi type. Ater substituting Askey- Wilson parameters given by ( [4. 14| ) and 
dividing by a constant with value (1 — g)^'' the rth difference equation in Theorem E 



for the Askey- Wilson type polynomials goes for g ^ 1 over in a differential equation 
^rPx ~ ^rxPx '^^ order 2r for the multivariable Jacobi type polynomials ||D1| , Sec. 
4]. The computation of the eigenvalue = limg^i(l — q)~'^'^E^Y hinges again on 
Lemma and the result is 



E^^ = Er {{p( + Al)^ ...,{pi + XnY; {pi)\ (piY) (5.20) 

with Er{- ■ ■ ; ■ ■ ■ ) taken from ( ^.6|) and 

p-; = {n-j)u + {iyo + iy,)/2, (5.21) 

where v = g ■, = go + q'q and vi = gi + g'l- For r = 1 the differential operator 

= limg^i(l — q)~'^^'D:^^ is given by D'^ ( p. 221) . More generally one has that 
is of the form 

^r= E (5.22) 

JC{1 n} jej 

\J\=r 

(where l.o. stands for the parts of lower order in the partials), but it seems difficult 
to obtain the relevant differential operators for arbitrary r in explicit form starting 
from D^^ fO). 



Theorem 5.5 (| pi| . Sec. 4]). The multivariable Jacobi polynomials p'l, A G A 
satisfy a system of differential equations of the form 

Dipi = E'^^^pl r = l,...,n, (5.23) 

where D;^ = limq^i(l — q)^'^^D^^ is of the form (|5.22| ) and the corresponding eigen- 
values are given by ( ^.20|) . 

For r = 1 the Differential equation (|5.23|) coincides with the second order Differ- 
ential equation (|3.24|) . 

Remarks: i. The proof in [pl| , |D4|] demonstrating that the multivariable Askey- 
Wilson type polynomials satisfy the system of difference equations in Theorem |5 
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runs along the same lines as the proof for the special case when r = 1 (Proposi- 
tion |3.1|) : it consists of demonstrating that the operator is triangular with 
respect to the monomial basis {Tn^^}x£A and that it is symmetric with respect to 
the inner product (■,-)a'4w- The triangularity proof consists of two parts. First it 
is shown that D 

AW^AW 

lies in the space Ti^^ by inferring that poles originating 
from the zeros in the denominators of v"^^{z) ( 15. 3| ) and w"^^{z) ( ^.4| ) all cancel each 
other. Next, it is verified that the operator indeed subtriangular by analyzing 
the asymptotics for x at infinity. 

It is noteworthy to observe that if one tries to apply the same approach to arrive 
at a direct proof of the system of difference equations for the multivariable Wilson 
and continuous Hahn type polynomials (i.e., without using the fact that these two 
cases may be seen as limiting cases of the Askey- Wilson type polynomials), then 
some complications arise. In both cases the proof that the operator Dr maps the 
space TC into itself and that it is symmetric with respect to the inner product (■, ■)a 
applies without significant changes. However, it is now not so simple to deduce from 
the asymptotics at infinity that the difference operator D,, is indeed triangular. The 
reason — in a nutshell — is that the functions w^{z) (|5.1CI|) and w'^{z) (|5.16| ) no longer 
have constant asymptotics for z — >• cxo (as does w"^^{z) for exp{iaz) oo). As a 
consequence, it is now much more difficult than for the Askey- Wilson case to rule out 
the a priori possibility that monomials with /i ^ A enter the expansion of D^nix 
in monomial symmetric functions (in principle Dr might a priori raise the degrees of 
the polynomials). For r = 1 one easily checks by inspection that such with ^ A 
indeed do not appear in the expansion of Drm\ but for general r this is not so easily 
seen from the explicit expressions at hand. 

a. Some years ago, a rather explicit characterization of a family of commuting 
differential operators simultaneously diagonalized by the Jacobi type polynomials 
(and generating the same algebra as the operators D(, . . . , D^) was presented by 



Debiard |De]. Alternatively, it also turned out possible to express such differential 



operators in terms of symmetric functions of Heckman's trigonometric generalization 



of the Dunkl differential- reflection operators related to the root system BCn |0 |Du |. 
In both cases, however, it is a nontrivial problem to deduce from those results an 
explicit combinatorial formula for (the coefficients of) the operator for arbitrary 
r. In the case of Debiard's operators one problem is that the corresponding eigen- 
values do not seem to be known precisely in closed form (and also that one would 
like to commute all coefficients to the left); in the case of an expression in terms of 
trigonometric Dunkl type differential-reflection operators it appears to be difficult to 
explicitly compute the differential operator corresponding to the restriction of the 
relevant symmetrized differential-reflection operators to the space of symmetric poly- 
nomials (except when the order of the symbol is small). In the latter case the problem 
is that one has to commute all reflection operators to the right (and preferably all 
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coefficients to the left). This poses a combinatorial exercise that seems tractable only 
for small order of the symbol. 

in. If one transforms the operator D'^ = D( to a second order differential operator 
that is self-adjoint in L^(R", dxi ■ ■ ■ dxn) by conjugating with the square root of the 
weight function A'^(a;) ( p.l9| ), then one arrives at a Schrodinger operator of the form 

(A-^)-3D^(A^)i = (5.24) 
^ 52 ^ -a^ {^^^^^ ~ ' ^^^^^ " 



J A ^-^ \ — i) COS"'' — ■■ 



2 



2 ' ' " i^f:nf^„^sin^ %{^3 - ^k) ' sin^ f (xj - Xk] 



i<j<k<n 

„2 



where = ct X]i<j<n(p/) • Explicit expressions for n independent commuting differ- 
ential operators generating the same commutative algebra as transformed operators 
r = 1, . . . ,n can be found in the literature as a special case of 
the formulas presented in ||00S| , |0S|| . However, just as in the previous remark it is 



again nontrivial to determine explicitly the relation between our transformed opera- 
tors (A-^)-2L)/(A-^)2, r = 1, . . . , n and the relevant specialization of the differential 
operators in [pOS| , |0S|| , because to our knowledge the eigenvalues of the latter oper- 
ators are not available in closed form (except in cases when the order of the symbol 
is small). 

6. Recurrence relations 

In this section we first recall the system of recurrence relations for the multivariable 
Askey- Wilson type polynomials presented in | |D5| ]. Next, the limit transitions of 



Section ^ are applied to arrive at similar systems of recurrence relations for the 
multivariable Wilson, continuous Hahn and Jacobi type polynomials, respectively. To 
describe these recurrence relations it is convenient to pass from the monic polynomials 
P\{x) to a different normalization by introducing 

Pxix)^c,p,{x), c, = c\'\^^±^ (6.1) 

A+(p + A) 

where c denotes some constant not depending on A (recall also that |A| = Ai + - ■ ■ + Xn) 
and the function A+(x) is of the form 

^+{x) = Y\. dv,+ ixj + Xk) dv,+ {xj - Xk) Y\. dw,+{xj)- (6.2) 

l<j<k<n l<i<" 

The precise value of the constant c, the vector p = (pi,... ,p„) and the form of 
the functions dv,+, dw,+ depends on the type of polynomials of interest and will be 
detailed below separately for each case. 
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The general structure of the recurrence relations for the renormalized polynomials 
P\{x) reads 

Er{x)Px{x)= Yl Ujc^,_^j\{p + \)VeJ,Ap+X)Px+e,j{x) 

.JC{1,... ,n},0<|J|<r (6-3) 

ej=±l, j&J; e^j+AeA 

with r = 1, . . . , n, and 

Ve.J,K{x) = Y[w{ejXj) Yl v{ejXj + SyXy) v{ejXj + Ej'Xj' + 1) 

j<j' 



Y\ v{6jXj + Xk) v{€jXj - Xk), 



UkA^) = 

(-l)P ^ (j^ w{eixi) Y\ v{eixi + evxv) v{-eixi - £i>xv - 1) 

ei=±l, leL Kl' 

X Y\_ v{eixi + Xk)v{eixi - Xk)y 

keK\L 

The functions -E'i(x), . . . ,En{x) appearing in the l.h.s. of ( |6.3| ) denote certain (ex- 
plicitly given) symmetric polynomials that generate the algebra of all symmetric 
polynomials. The rth recurrence relation expresses the fact that the expansion of the 
product Er{x)Px{x) in terms of the basis elements P^^), /i G A is known exphcitly 
(i.e., the coefficients in the expansion are known in closed form). The expansion 
coefficients are determined by the functions v and w (together with the vector p) 
whose precise form again depends on the class of polynomials of interest. In com- 
binatorics one sometimes refers to this type of recurrence relations as to generalized 
Fieri type formulas after similar expansion formulas for the products of elementary or 
complete symmetric functions and Schur functions (in terms of the latter functions) 
[M4|| . In the simplest case, i.e. when r = 1, the Recurrence formula ( |6.3| ) reduces to 
an expression of the form 

Eix)Pxix) = Yl V,ip + X){Px+,^ix)-Pxix))+ (6.4) 

l<j<n 



l<j<n 
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with E{x) = Ei{x) and 

V±j{x) = i^i^Xj) Y\_ ^i^^j + ^k) v{±Xj - Xk)- (6.5) 

l<fc<n, k^j 

We will see that for the four families considered below this formula reduces in the case 
of one variable to the well-known three-term recurrence relation for the Askey- Wilson, 
Wilson, continuous Hahn or Jacobi polynomials, respectively. 

Note: To date a complete proof (contained in |P5 |) of the recurrence relations for 
the multivariable Askey- Wilson type polynomials exists only for parameters satisfying 
the (self-duahty) condition 

qtotiW = l. (6.6) 

It is for this reason that in the present formulation of the system of recurrence rela- 
tions for our multivariable hypergeometric polynomials given below some restrictions 
on the parameters are imposed (except in the case of Jacobi type polynomials where 
such a restriction turns out to be not necessary). As is argued in [ P5|| for the Askey- 
Wilson type, however, it is expected (and easily checked for the special case of one 
variable) that our recurrence relations remain valid also for parameter values not 
meeting these restrictions (cf. Remark i. at the end of this section). 

6.1. Askey- Wilson type. The normalization constants for the polynomials in 
( |6.1| ) of Askey- Wilson type P^^{x) are determined by the constant = 1 and the 
functions 

<^(.) = t-^^'^^, (6.7) 

[t^q^, tig^ t2g^, hq""] q)oo 

Here we have introduced dependent parameters ir that are related to the Askey- 
Wilson parameters tj. by 

to = {W2hq-^fl\ 

ii = {totit2^t^^qy/'^, , . 

is = {totiH^%qy/^ 

The vector p = p'^^ = {pf^^ , ■ ■ ■ , Pn^) has in the present case the components 

pj^ = Mogr,, r, = e-Kt,t,t2hq-'f'^ (6.10) 

and is introduced mostly for notational convenience. In fact, the logarithm with base 
q entering through the components pf^ merely has a formal meaning and appears 
in our formulas always as an exponent of q. 
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The symmetric functions Ef^{x) multiplying P^^{x) in the l.h.s. of (|6.3| ) are 
given by 

Ef'^ix) = Erie''''' + e-*°"S . . . , e'"^'" + e"'"""; + f-\ . . . , f„ + f"^) 

(6.11) 

with Er{- ■ ■ ; ■ ■ ■ ) being taken from (|5.6|) (as usual) and with 

f, =r-^(totitW)'/' (6.12) 

(cf. ( 6.10 )). The coefficients in the r.h.s. of ( |6.3|) entering the expansion of the 
function E:^^ (x) (x) in terms of the basis elements P^^{x) are characterized by 
the functions 

v^^iz) = (6.13) 

-AWf \ tl Z Z 1 -l\-l/2 no<r-<3(-'- ~ ^'^ 5'^) ( a -\ A\ 

w (z) = {W^hq ) ' (^_-^2;)(i_g2.+i) (6-14) 



(cf. ( |5.3| ), (|5.4|) ). We now have the following theorem from |P5| 



Theorem 6.1 ( [p5|| ). The renormalized multivariable Askey-Wilson type polynomi- 
als P^^{x), A G A (|2.2| ), satisfy a system of recurrence relations given by (|6.3|) (with 
Ef^, v^^, w^^ and p^^ taken from ( |6TID , ( lOl , ( lOl and { WM , and with 
r = 1, . . . , for Askey-Wilson parameters subject to the condition qto t]~^'^2 ^^3 ^ — ^■ 

For r = 1 one has 

E^^{x) = E^^{x) = J2 (e'""' + - ^1 - ^,"') • 

i<j<" 

The corresponding Recurrence formula (|6.4|) , (|6.5| ) coincides in the case of one variable 
with the three-term recurrence relation for the renormalized Askey-Wilson polyno- 
mials 

prw ^.^3 ('"•'°«;f;;'f7-'°^"" (6.15) 

V tO^l, t^ol2, tQlS / 



which reads p| 

(2 cos(ax) - to - to ^)^i^^(a;) = (6-16) 
(1 - totit^t^q'-') ni<r.<3(l - toW) 



toil - totit2tsq^'-^)il - totit2tsq^' 
toil-q')Ui<r<s<3i^-isW'' 



(p,^f (x) - pn^)) + 



(1 - totit2t3q^'-^){l ~ totit^tsq^^- J 
It is clear that in this special situation Theorem |6.1| indeed holds without restriction 



on the parameters as was anticipated more generally for arbitrary r and n in [D5 
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It is important to observe that, despite the appearances of square roots (and infinite 
products) in intermediate expressions, both the normahzation constants ^ and the 
Askey- Wilson type Recurrence relations ( |6.3|) are in the end rational in q and the 
parameters t, to, . . . , ts- For the recurrence relations this is is rather immediate from 
the fact that 



V 



^^'e,{pr + ^i) + ^dpr + A.)) = t-"^ I I (6.17) 

-AW/ / AW I \ \\ /i i 1 i -lN-1/2 no<r<3(-'- ~ tr^ ^ ^) /r, 

^ +\)) = {W,t,q ) (^_^,-^-,.^..^(^_^,.^^,.^.^._,i) (6.18) 

The point is that the combinations t^W^'', tr t^^ and r^^.? are rational in t, t^, . . . , ^3 
plus that the functions always emerge in pairs in ( |6.ij| ) and (totit2^3?~^)^^^ = ^0 
(hence the square root constant factors in , do not spoil the rationality). The 
last equality is also needed to see that E:^^{x) ( |6.11| ) in the l.h.s. of the recurrence 
relation depends rationally on the parameters: fj = t"'~-' (to^i^2^3Q''^)^^^ = t^~-'to- For 
the rationality in the parameters is seen similarly after rewriting in the form (by 
cancelling common factors in numerator and denominator) 

i<i<fc<n ^)^^+^^- (l)x,~\k 



AW 



n 



l<i<n ^^O^i' ^l^J' ^STj; q)xj 



where C = Ui<j<n^j 



6.2. Wilson type. If we substitute Askey-Wilson parameters in accordance with 
( [4.7| ), then the polynomials P^^{x) converge for a — > to renormalized multivariable 
Wilson type polynomials P^{x) = cYpx{x). The normalization constant cY is of the 
form in ( |6.1|) , (|6.2|) with = —1, the vector taken from ( 5.12 ), and with the 
functions d^_^_, given by 

= = (6.20) 

Here we have introduced dependent parameters Ur related to the Wilson parameters 
Ur by 

i^O = (t'O + 1^1 + 1^2 + 1^3- l)/2, 

z>i = (z/o + z/i-z/2-z^3 + l)/2, .g2i) 

Z>2 = (Z/Q - '^l + ^^2 - «^3 + l)/2, 
^^3 = (t'o - '^l - ^^2 + ^^3 + l)/2. 
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To verify that the Askey- Wilson type polynomials P^^{x) with parameters ([4.7| ) 
indeed converge for a — ^ to the renormalized Wilson type polynomials P'^{x) thus 
defined, one uses Proposition and the fact that (for parameters ([4.7|)) 

lim(-a2)l^lc^^ = cf. 

The latter limit can be easily checked with the aid of Representation (|6.19|) for 
entailing 

( W \ W\ I W W\ 



rv — >(] 

(2pr 



l<j<n + ' ^1 + ' + ' ^3 + Pj )X, 

It is not difficult to see that the r.h.s. of this expression indeed coincides with the 
above defined cY of the form (|6.1| ), ( |6.2|) by rewriting all Pochhammer symbols as a 
quotient of two gamma functions. 

Let us next turn to the limiting behavior of the recurrence relations. For parameters 
as in (4/7) and after division by a"^^ the rth recurrence relation for the Askey- Wilson 
type polynomial P^^{x) goes over in a recurrence relation of the form in (|6.3| ) for 
the Wilson type polynomial P^{x). The relevant symmetric function in the l.h.s. of 
the resulting Wilson type recurrence relation is given by 

4^(^) = (-ir E n^? E (pY---pY-J' (6-22) 



,/c{i,... ,n} j£j r<li<---<lj._i ji<n 
0<\J\<r 



with 



Pj =in- j)iy + (i)o + z/i + z/2 + z>3 - l)/2 (6.23) 



(cf. ( ^.12|) ). Furthermore, the coefficients in the r.h.s. of this recurrence relation are 
determined by the vector (|5.12| ) and the functions 

z 2z{2z + 1) 

(cf. (|5.1CI|) ). To check this transition from the Askey- Wilson to the Wilson type 
recurrence relations one uses that for Askey- Wilson parameters given by (|4.7] ) one 
has 

lim a-^'Ef^ix) = EY{x) 
(in view of Lemma b.2\) and that 



lim^ v^^ie.pf^ + SkPr + ^)= v'^ie.pj + e^p^ + z), 



1- -2 ^AW/ AW I \ ^W/ W , \ 

lim a w £iPf + z) = w (€iPY + z) 
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We have thus derived the following theorem. 

Theorem 6.2. The renormalized multivariahle Wilson type polynomials [x), A G 
A (pl^), satisfy a system of recurrence relations given by ( |6.3|) (with , -0^, 
and taken from ( |6.22| ), (|6.24]) and ( ^.12|) , and with r = 1,... ,n) for Wilson 
parameters subject to the condition Uq — Ui — 1/2 — 1^3 + 1 = 0. 



The condition on the parameters in Theorem |6.2| is an immediate consequence of 
the parameter condition in Theorem and may be omitted once we are sure that 
Theorem |6.1| holds for general parameters. 

For r = 1 the recurrence formula is given by (|6.4D, (|6.5|) with 



E^ix) = E^ix) = -Y: {^] + {Pjf) ■ (6.25) 

l<i<n 

In the case of one variable this recurrence formula coincides with the three-term 
recurrence relation for the renormalized Wilson polynomials 

j^w, X f-^^ uo + ui + U2 + U3 + I - 1, uo + ix, uo - ix \ 

F, (x) = 4F3 ; 1 ) 

V Uo + Ul, + 1^2, 1^0 + ^3 ) (g 26) 

which reads |p<.l5|| 

- (x^ + vl)Pr{x) = (6.27) 

{I + Uq + Ui + U2 + - 1) ni<r<3(^ + 1^0 + l^r) 



{pr^,{x)-pr{x)) + 



{21 + UQ + U1 + U2 + U3 - l){2l + PQ + U1 + U2 + U3) 

^ ni<.<.<3(^ + ^s + Vr - 1) _ . 

{21 + vo + vi + V2 + v-i- 2){2l + UO + UI + 1^2 + 1^3-1)^ I y !)■ 

This three-term recurrence formula indeed holds without restriction on the parame- 
ters. 



6.3. Continuous Hahn type 

( [4.11| ) and shift the variables 



X 



the Askey- Wilson type polynomials P 
tinuous Hahn type polynomials P 



If we substitute Askey- Wilson parameters given by 
,■ over a half period (cf. ( [4.10D ), then for a -h> 



AW, 



X 



cHi 
A 



c'^^p'^^{x) of the form in (|6 



the vector p'^^ taken from ( |5.18 



become renormalized multivariable con- 

, (13) with 
d'^, reading 



and with the functions d^^j^, 



T{iy + z) 

r{z) 



n 



0<r<2 



r(z>, + Z) 



T{2z) 



(6.28) 
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Here we have again employed dependent the parameters Oq, . . . ,1)3 that are now 
related to the Wilson parameters z/^, z/f by 

= {1^0 + i^t + 1^0 + 1^1 - ^)/'^, 

z>2 = - ^1+ - + + 1)/2, ^^•^''^ 

(Notice, however, that df^j^ depends only on z>o, z>i, z>2 and not on z>3, which is merely 
introduced for convenience and will be used below (cf. ( |6.31|) ).) To verify the tran- 
sition P^^ (x) P^^{x) one uses Proposition [4.2| and the fact that (for parameters 
( Pl|) ) lim„_o(2a)l^lcf'^ = cf . The derivation of the limit (2a)l^lcf^ is very 

similar to that of the Wilson case and hinges again on Representation (|6.19|) : first 



one writes cf^ ( 6.19| ) explicitly as a rational expression in the Askey- Wilson param- 



eters t, to, . . . ,^3 by invoking of the definitions ( |6.9|) and ( |6.10|) ; next, substituting 
the parameters ( [4 .111 ) and sending a to zero after having multiplied by (2q;)''^' entails 



an expression for cf involving Pochhammer symbols 



l<j<k<n + Pf + )a.+A, (^ + Pf - Pf)x,-X, 



X 



n 



cH\ 



l<,<n (^0+pf, ^2+pf)A, 



which is seen to be equal to the stated expression of the form (|6.1|), ( |6.2|) by rewriting 
all Pochhammer symbols as a quotient of gamma functions. 

The corresponding recurrence relations for the renormalized multivariable Wilson 
type polynomials P^^{x) are of the form given by ( |6.3|) with the symmetric functions 
in the l.h.s. given by 

Ef{x) = {-IY Yl U'^^ E (6-30) 



JC{1,... ,n) j£j ^■<h<---<'r-|J| 

0<\J\<r 



with 



pf = in- j)iy + (i)o + z>i + z>2 + ^>3 - l)/2. (6.31) 



The coefficients in the r.h.s. are determined by the vector p'^^ ( ^.18|) and the functions 

V {z) = , w [z) = -- ■ 6.32 

z 2z[2z + 1) 

To arrive at these recurrence relations for the multivariable continuous Hahn type 
polynomials we have substituted the parameters ( [4.11|) , shifted the variables over a 



half period ( [4.10|) and divided both sides by (2m)'". For a — the rth Askey- Wilson 
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tj^e recurrence relation then passes over in the rth continuous Hahn type recurrence 
relation. At this point one uses that (for parameters as in ( |4.11| )) 

hm (2m)-'^4^^(x-f ))=4^^(x) 

(which immediate from the definition of E^"^ [x)) and that 

lim i)^^(e,p/^ + e.pr + ^) = v'^''{e,pf + e,pf + z\ 

lim {2iaY^w^'^{e,pf ^ z) = w'^'ie.pf + z). 

We arrive at the following theorem. 

Theorem 6.3. The renormalized multivariable continuous Hahn type polynomials 
P^^{x), A G A (p.2|), satisfy a system of recurrence relations given by ( |6.3| ) (with 
i^cH ^ -cH ^ ^cH pcH ^Q/jgj^ from ( |6.30D , ( |6.32| ) and ( |5.18| ), and with r = 1, . . . ,n) 
for continuous Hahn parameters subject to the condition Vq — — Vq — + 1 = 

The condition on the continuous Hahn parameters VQivf in Theorem p.3| is of 
course an artifact of the parameter condition in Theorem |6.1| and may again be 
omitted as soon as it is shown that Theorem ^Tl] holds for general parameters. 

For r = 1 the recurrence relation is of the form in (16.41), (B.SI) with 



E^"{x) = Ef{x) = - [ixj + pf) . (6.33) 

i<j<" 

In the case of one variable this formula coincides with the three-term recurrence 
relation for the renormalized continuous Hahn polynomials 



r^cHi N _ , -I, + + + Vi + I - I, + ix 



(6.34) 



which reads 



-{ix + u^)Pr{x) = (6.35) 

(/ + U+ + Z/+ + Z/q- + Z/f - !)(/ + Z/o+ + + + Z/f ) ( .H,. _ pcH(^\ , 

[21 + + u{ + + - l)[2l + + u{ + + v^) ^ 

-Ki + + - + + - 1) (p-^ix) - p'^'ixW 

{21 + v+ + Z/+ + z/q- + z/^- - 2)(2/ + z/o+ + Z/+ + z/q- + z/f - 1) ^ '"^^ ^ ' ^ 
and holds for general parameters. 

6.4. Jacobi type. After substituting the Askey- Wilson parameters ( [4.14|) and send- 
ing q to one, the Askey- Wilson type polynomials P^"^ {x) pass over to renormalized 
multivariable Jacobi type polynomials Px{x) = c'l p'l{x) of the form in ( |6.1D , (|6.2|) 
with the Jacobi parameters and z/i taking the value Qo + Qq and gi+g'i, respectively. 
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The normalization constant cl is determined by the constant c"^ = 2 ^ and the vector 
p"^ taken from ( |5.21| ) together with the functions d'^ and (i;^ _,_ given by 



T(z) 



di.Az) 



V{u^ + z)V{Pi + z) 
Vi2z) 



(6.36) 



Here we have introduced dependent parameters z>o, z>i related to the Jacobi parame- 
ters 1^0! by 



z>i = (z/o - z/i + l)/2. 



(6.37) 



The verification of the transition P^"^ [x) 



— * P\{x) hinges on Proposition [4.3| and the 
fact that (for parameters ( |4.14| )) limg^i = cj^. The latter limit is again checked 
using Formula ( |6.19|) for cj^^ entailing for g — 1 



2 2A| JJ 



{pi - pi] 



At \h 



l<j<k<n 



X 



n 

l<j <n 



(z/ + p/ + pi)x^+x, {y + p/ - pi)\,^\^ 

(2p/)2A, 



(z>o + p/, z>i + p/)a, 



which seen to be equal to the stated expression for of the form ( |6.1| ), ( |6.2| ) by once 
again rewriting the Pochhammer symbols in terms of gamma functions. 

The corresponding recurrence relations for P'l{x) are of the form in ( |6.3D , with the 
symmetric function in the l.h.s. given by 



Ei(x) 



-1^ 



JC{1,... ,n) jeJ 
\J\=r 



(6.38) 



(this is (—1)*^ times the rth elementary symmetric function in the variables sin {-^), 
j = 1, . . . ,n) and the coefficients in the r.h.s. being determined by p'^ ( [5.21[ ) and the 
functions 



v-^{z) 



V + z 



z 



w\z) 



[uq + z){i)i + z) 
2z{2z + l) 



(6.39) 



One obtains the recurrence relations for the multivariable Jacobi type polynomials 
for q —>■ 1 from the Askey- Wilson type recurrence relations with parameters ( |4.14| ) if 
one divides both sides by 2^^'. To check this one uses that (for parameters given by 

(EB) 



2r ibAW , 



hm 2-'''E. 

q-^l 



EHx) 
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(which is rather immediate from the definition of Ef^{x)) and that 
hm {)^^(£,p/^ + Skpt"^ + z)= v-\e,p'l + eupi + z) 

hm 2-^w^^{e,p^^ + z)= w\e,p-\ + z). 

We thus arrive at the foHowing theorem. 



Theorem 6.4. The renormalized multivariahle Jacobi type polynomials P/(x), A G 
A ( |2.2| ), satisfy a system of recurrence relations given by ( |6.3| ) (with , ■0"', w"^ and 



p^ taken from (|6.38|) , ( |6.39|) and ( |5.21|) , and with r = 1, . . . ,nj. 



Notice that in the Jacobi case it is not needed to impose any condition on the 
parameters z/q, vi. The point is that if we substitute the Askey- Wilson parameters 
( f4.14| ) then the condition q to t]~^^2 ^ = 1 in Theorem ^Tl] gives rise to the condition 



9o — 9'q — 9i~ 9i = on the parameters go, gi, g^, g[. However, given this condition the 
Jacobi parameters i^o = go+g'o and = gi+g'i can still take arbitrary values. In other 
words: the confiuence of the parameters in the limit g ^ 1 has as consequence that the 
full three-parameter family of multivariahle Jacobi type polynomials (parametrized 
by z/, z/q and z/i) may be seen as limiting case of the family of (self-dual) Askey- Wilson 
type polynomials with parameters satisfying (|6.6|) . 

For r = 1 we now have a recurrence formula of the form in (16.41) with 



E-^ix) = E({x) = -J2 si^' (^) • (6-40) 

In the case of one variable this recurrence formula reduces to the three-term recurrence 
relation for the renormalized Jacobi polynomials 

,J/ N r f~^^ Uo + l^l + l _ . 2f(^^\ 



reading ||AS| , |Kl5|| (notice, however, that our normalization of the polynomials differs 
slightly from the standard normalization) 

-sin2(^) P,\x)= (6.42) 

(/ + z/o + z.O(/ + z.o + l/2) _ 
(2/ + z.o + ^i)(2/ + z.o + z/i + l) ^^'^^^"^^ ^'^"^^^^ 

- ^1±^:LzV^) (P^(x) -P/(x)) . 

Remarks: i. The combinatorial structure of the recurrence relations for the mul- 
tivariahle Askey- Wilson type polynomials is very similar to that of the difference 
equations in Section |^. This is by no means a coincidence. In fact, in [ p5|] the 



recurrence relations were derived from the difference equations with the aid of a 
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duality property for the renormalized multivariable Askey- Wilson polynomials first 
conjectured by Macdonald 

Here P^^{x), /i G A, denotes the renormalized multivariable Askey- Wilson type 
polynomial with the parameters to, ■ ■ ■ , ^3 being replaced by the parameters to, . . . , ^3 
(cf. (|6.9|) ) and p"^^ is given by ( |6.10| ) whereas is the corresponding vector 



with Tj ( |6.10D replaced by fj ( |6.12| ). If one substitutes x = ^ (p + A) in the 
rth difference equation of Theorem |5.1| for the polynomial P^'^(x), then by using 



Property ( 6.43| ) (and the fact that the coefficients have a zero atx = ^(p^^ + A) 



if A + Cej 7^ A (|2.2|) ) one arrives at the rth recurrence relation of Theorem |6.1| (first 
at the points x = ii£M(p^'^ _)_ ^ ^ and then for arbitrary x using the fact that 
one deals with an equality between trigonometric polynomials). In [ p5|| we proved 



Property ( |6.43D (and hence the recurrence relations) for parameters satisfying the 
self-duality condition (|6.6|) (implying that ir = tr)- It is of course expected, however, 
that Macdonald's conjecture (|6.43| ) holds for general parameters. In the case of one 
variable ( |6.43D this is immediate from the explicit expression of the polynomials in 



terms of the terminating basic hypergeometric series ( |6.15 ) 



a. By applying the transition Askey- Wilson — > Wilson to (|6.43| ) one arrives at a 



similar duality relation for the multivariable Wilson type variables 

Pf (z(p^ + p)) = (z(p^ + A)) , (6.44) 

where P^(a;) now denotes the Wilson type polynomial with the parameters z/q, . . . ,1^3 
being replaced by z>o, . . . , z/3 (cf. ( |6.21| )). The self-duality condition for the parameters 
inherited from ( |6.6| ) then becomes z/q — z/i — z/2 — ^3 + 1 = (cf. Theorem |6.2| ), which 
implies that Vr = i^r- For n = 1 the relation ( |6.44| ) is again immediate for general 
parameters from the explicit expression of the polynomials in terms of the terminating 
hypergeometric series ( |6.34D . 

in. For p = the r.h.s. of ( |6.43D becomes identical to one and one arrives at 
the relation P^"^ {ia'^^hif) = 1 (where Inf stands for the vector (Infi,... ,lnf„) 
with fj taken from (|6.12|) ). This equality amounts to the following evaluation or 



specialization formula for the monic Askey- Wilson type polynomials (cf. Definition 
□)) 

[la mr = — ^7 . 6.45 
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Using the limit transitions of Section ^ one arrives at similar evaluation formulas for 
the multivariable Wilson, continuous Hahn and Jacobi type polynomials: 



p^(o) 



) = (-i)i'' 



A!^(p^ + A) 



22|A| 



Aiip-^ + X) 



Hip') 



(6.46) 
(6.47) 
(6.48) 



For the Askey- Wilson case this evaluation formula was conjectured by Macdonald 
and then proven by Cherednik [|C2|| for special parameters related to reduced root 
systems (and admissible pairs of the form {R,R}^)) with the aid of shift operators. 
A similar proof for the Jacobi case can be found in ||0p|| . The proof of the Askey- 
Wilson type specialization formula ( |6.45| ) based on the recurrence relations for the 
more general parameters subject to the condition (|6.6D was first presented in ||D5|. 



7. Orthogonality and normalization 

The difference/differential equations in Section ^ hold for generic complex param- 
eter values in view of the rational dependence on the parameters (cf. also Remark ii. 
of Section The same is true for the recurrence relations in Section ^ (although we 
still have to impose the parameter restrictions stated in the theorems of Section ^ so as 
to keep our present derivation of the recurrence relations to remain valid, cf. Remark 
i of Section |^). In order to interpret the polynomials as an orthogonal system with 
weight function A, however, rather than allowing generic complex parameters we will 
from now on always choose the parameters from the domains given in Section 0. (The 
conditions on the parameters in Section ^ ensure that the relevant weight functions 
A{x) are positive and that the integrals defining the associated inner products (■, ■)a 
converge in absolute value.) 

It is immediate from the definition in Section ^ that the polynomial px is orthogonal 
to for /i < A. It turns out that also for weights that are not comparable with respect 
to the partial order ( ^^.3] ) the associated polynomials are orthogonal. 

Theorem 7.1 (Orthogonality). The multivariable Askey- Wilson, Wilson, continu- 
ous Hahn and Jacobi type polynomials p\, A e A, form an orthogonal system with 
respect to the inner product with weight function A{x), i.e. 

(Pa,Pm)a = if X^p (7.1) 

(for parameters with values in the domains given in Section\^. 

The proof of this theorem hinges on the following proposition. 
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Proposition 7.2 (Symmetry). The operators Di, . . . , of Section ^ are symmet- 
ric with respect to the inner product (-, ■)a, i-c. 

{Drmx,m^)A = {mx,Drm^)A (7.2) 

(for parameters with values in the domains given in Section^. 

A detailed proof of the symmetry for the Askey- Wilson type difference operators 
can be found in [pl| , Sec. 3.4]. The same reasoning used there can also be applied 



to prove the proposition for the Wilson and continuous Hahn case. The Jacobi case 
follows from the Askey- Wilson case with the aid of the limit transition from Askey- 
Wilson type to Jacobi type polynomials and operators (we refer to [ [Pl| , Sec. 4] for 
the precise details). 



By combining Proposition |7]2| with the results of Section ^ we conclude that 



the polynomials p\, A G A, are joint eigenf unctions of n independent operators 
Di, . . . ,Dn that are symmetric with respect to (■, ■)a- The corresponding (real) 
eigenvalues Ei^x, . . . ,En,x separate the points of the integral cone A ( ^^ , i.e., if 
Er^x = -E'r.At for r = 1, . . . ,n then A must be equal to fi (this is seen using the fact 
that the functions Er{xi, . . . , x^, Vr, ■ ■ ■ Un) (|5.6| ) generate the algebra of permutation 
symmetric polynomials in the variables Xi, . . . ,a;„). It thus follows that the polyno- 
mials px and with A 7^ /i are orthogonal with respect to the inner product (■, ■) a as 
eigenfunctions of a symmetric operator corresponding to different eigenvalues. An- 
other proof for the orthogonality of the multivariable Askey- Wilson type polynomials 
was given by Koornwinder in 0. In the case of Jacobi type polynomials independent 
orthogonality proofs (not using the fact that the Jacobi type is a degenerate case of 
the Askey- Wilson type) can be found in e.g. [De] or in Q (upon specialization to 
the root system BCn)- 

The main purpose of this section is to express the squared norms of the polynomials 
(viz. {px,P\)a) in terms of the squared norm of the unit polynomial (viz. (1, 1)a, 
which corresponds to A = 0). Our main tool to achieve this goal consists of the 
recurrence relations derived in the preceding section. Starting point is the identity 

{ErPx, i'A+a;,)A = (^A, A+u;.)a, 0;^ = Ci H h 6^, (7.3) 

where Px{x) denotes the renormalized Askey- Wilson, Wilson, continuous Hahn or 
Jacobi type polynomial of the form ( |6.1| ) and Er{x) is the corresponding symmetric 
function multiplying Px{x) in the l.h.s. of the Recurrence relation ( |6.3| ). (In all four 
cases of interest the functions Er{x) are real for parameters in the domains of Section]^ 
and we hence have ( |7.3| ) trivially.) If we work out both sides of (|7.3|) by replacing 
ErPx and E^Px+uir by the corresponding r.h.s. of ( |6.3|) and use the orthogonality 
of the polynomials (Theorem |7. 1|) , then we arrive at the following relation between 
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{P\,P\)a and {Px+u,,.,P\+iVr)A 

V{1,... ,r},{r+l,... ,n}ip + >^){Px+u,r, Px+iVr) A = (7.4) 

^{l,...,r},{r+l,...,n}(— P — A — UJr){Px, Px)a- 

(Recall that Uk,p = 1 for p = and that V{i....,r}.{r+i,...,n} is taken to be in accordance 
with the definition in ( |6.3|) with all signs ej, j & J being positive.) In principle one 
can use this relation to obtain for each A G A the squared norm of P\{x) in terms 
of the squared norm of Pq{x){= 1) by writing A as a positive linear combination 
of the (fundamental weight) vectors tui, . . . ,ujn and then apply (|7.4| ) iteratively by 
walking to the weight A starting from the zero weight (0, . . . ,0) through the successive 
addition of fundamental weight vectors uJr- Indeed, it is not difficult to verify that 
the combinatorial structure of the coefficients V{i^ ,,^r},{r+i,...,n} is such that the result 
does not depend on the order in which the fundamental weight vectors uJr are added, 
i.e., the result is independent of the chosen path from (0, . . . ,0) to A. This hinges on 
the easily inferred (combinatorial) identity 

^{l,...,r},{r+l,...,n}(— a; — ^t) ^{1,... ,s},{<i+l,... ,n} (" 3^ " ^s) _ 

V{l,...,r},{r+l,...,n}(a; + UJs) ,s},{s+l,... ,n} (a^) 

V{l,...,s},{s+l,...,n}(— a; — UJr— ^^s) ^{1,... ,r},{r+l,... ,n} ( — 3^ — ^^r) 
y{l,...,s],{s+l,...,n]{x + iOr) V{i^,„ ^r},{r+l,... ,n}{x) 

which expresses the fact that the result for the quotient of (Pa+lj,.+ws7 P\+ujr+uis) a and 
(Pa) P\)a computed via ( |7.4| ) does not depend on the order in which ujr and Ug are 
added (as it clearly should not). 

To write down the answer for (Pa, Pa)a/ (1, 1)a for general A G A we introduce the 
functions 

A±(x) = J]^ dv,±{xj + Xk)d^,±{xj - Xk) W dw,±{xj), (7.5) 

with dv^±{z) and dw,±{z) (7^ 0) satisfying the difference equations 

dy^+{z + 1) = v{z) d^^+{z), d^_{z+l)=v{-z-l)d^_{z), (7.6) 
d^^+{z + 1) = w{z) d^^+{z), d^_{z + 1) = w{~z - l)dw^^{z), (7.7) 

where v{z) and w{z) are taken to be the same as in Section It is immediate from 
the difference equations (|7.6|) , (|7.7| ) that 

A+(x + UJr) 



A+(x) 

A_{x + UJr) 

A_(x) 



^{i,...,r},{r+i,...,n}(a;), (7.8) 

V{l,...,r},{r+l,...,n}(-a; - UJr)- (7.9) 



HYPERGEOMETRIC POLYNOMIALS IN SEVERAL VARIABLES 



37 



With the aid of the properties 
form 



8f) and (17.91) we can rewrite Relation ([7.41) in the 



{Px,Px)a 



A+(p+A) 
A_(p + A) 



{Px+LUr^ Px+OJr)^. 



(7.10) 



By using the fact that the fundamental weight vectors ui, . . . ,u!n positively gen- 
erate the integral cone A (^.21), one deduces from this equation that the quotient 



(Pa, -Pa) a A+(p + A)/A_(p + A) in the l.h.s. of ( [7.10|) does not depend on the choice 
of A G A. Comparing with its evaluation in A = then entails 



{Px,Px)a A_(p + A) A+(p) 



(1,1). 



A+(p + A) A_(p) 



(7.11) 



As we will see below (and is suggested by the notation), it turns out that the function 
A+(a;) in the present section coincides with that of Section |[ So, by combining (|7.11|) 
with ( |6.1|) we obtain 



{Px,Px)a 
(1,1)a 



|-2|A| 



A+(p + A)A_(p + A) 
A+(p)A_(p) 



(7.12) 



We will now list for each of our four families AW, W, cH and J the associated 
functions d^^±{z) and dw,±{z), and formulate the corresponding evaluation theorem 
for the quotient of {px,Px)a and (1, 1)a- The proof of this theorem boils in each case 
down to verifying that dv±{z) and dw ±{z) indeed satisfy the difference equations 
(PI) and (O)- 

In the case of multivariable Jacobi type polynomials the value of the integral for 
the squared norms {px,Px)a was computed (in essence) by Opdam |Pp|| with the 
aid of shift operators (see also [[HS|| ) . Formulas for the squared norms of the Askey- 



Wilson type polynomials were first conjectured by Macdonald ||M2 , CT3 ] and then 
proven by Cherednik [ pi|| for special parameters related to the reduced root systems 
(and admissible pairs of the form {R, R"^)) using a generalization of the shift operator 
approach. Recently, Macdonald announced a further extension of these methods to 
the case of general Askey- Wilson parameters |[M5| . 

In ||D5|| the present author combined a proof of ( |7. 12| ) using the recurrence relations 
along the lines sketched above with Gustafson's constant term formula [CJu, Ka\ for 
(1, 1)aaw, which led to an alternative derivation of Macdonald's formula for the value 



of (pr^ 



Px^)a^^ (at present in the case of parameters satisfying Condition (|6.6D) 



The (Askey- Wilson) case is included here just for the sake of completeness. Also for 
the Jacobi type polynomials our proof of ( [7.12|) with the aid of recurrence relations 
rather than shift operators provides, when combined with the previously derived 
expression for the constant term (1,1) 



Ml 



I , an alternative way to demonstrate 
the validity of the known evaluation formula for {px,Px)aj- 
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7.1. Askey- Wilson type. In the case of Askey- Wilson type polynomials the func- 

5 "'W,+ \ 



tions A-|-(a;) = A^{x) of the form in (|7.5|) are characterized by d^^{z) ( |6.71) , d^^[z 
(U) and 

MW( \ _ fz/2 (^^"^^9)00 (j.o\ 

dT-iz) = m2hq-y' \^_, . , If tVh f-i , (7-14) 

(to ti ^2 <l^\ h 1^ , 1)00 

(with ir given by (|6.9|)). Formula ( |7.12|) leads us to the following theorem for this 
case. 



Theorem 7.3 ([D5|). Let us assume Askey-Wilson parameters with values taken 
from the domain indicated in Section \2.1\ such that the recurrence relations of Theo- 
rem \6. 1\ hold. Then one has 

{pf^^pf^)^^^ _ Kf^{pAw + A) K^W{pA^ + A) 
(with given by (|6.1CI|) / 



(7.15) 



To complete the proof of the theorem it suffices to infer that the functions d^^[z) 

and d^^{z) indeed satisfy the corresponding difference equations of the form in (|7.6|) 
and ( |7. 7| ) ; a fact not difficult to deduce from the standard relation for the g-shifted 
factorials (a; g)oo = (1 — c^) (o-Q'; (l)oo- 

For n = 1 the r.h.s. of ( |7.15| ) reduces to an expression of the form 

{Pl iPl }aAw/{Pq ,Pq }aAW 



with (cf. ||AW| , [KS|| and recall our normalization in ( p.8| ) 

Ana'^ {totit2t3q 

(totiMsg'"^ q)i (g'+^; g)oo Uo<r<s<3itrtsq^; g)oo ' (y^^g) 



/^AW ^AW\ 47ra ^ {totitit^q^^-q). 



7.2. Wilson type. The relevant functions A±(x) = A^(x) of the form in ( |7.5| ) are 
determined by d^^{z), d^_^_{z) ( |6.2CI| ) and 

^ _ Tj-U + Z + l) _ no<.<3r(-^r + ^+l) 

r(^ + i) ' <^^oA^)- r(2z + i) 

(with Or given by ( |6.21D ). The difference equations of the form in (|7.6|) and (|7.7|) 
are easily verified using the standard functional equation r(;z + 1) = 2;r(2;) for the 
gamma function. The normalization theorem becomes in the present situation. 



(7.18) 
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Theorem 7.4. Let us assume Wilson parameters with values taken from the domain 
indicated in Section |^.^ such that the recurrence relations of Theorem \6.2{ hold. Then 
one has 

{pY,pY)a^ _ A^(p^ + A) A^jp^ + A) 
{1,1) AW Af(p^)A^(p^) 

(with given by ( |5.12 ) / 

For n = 1 the r.h.s. of (|7.18| ) reduces to 

/ W W\ II W W\ 
{Pl ,Pl )aw/{Po ,Po )aw 

with (cf. |P:<.S|| and recaU our normahzation in ( |2.13| )) 

W 4^ ^' no<r<.<3 ^i^r + l^s + I) 



(pr.pDAw 



{i^o + 1^1 + 1^2 + 1^3 + I - 1); r(z/o + Z/i + Z/2 + Z/3 + 21) ' i^^^^Q^ 



7.3. Continuous Hahn type. The functions A-|-(x) = A^(a;) of the form in ( |7.5| ) 
are now determined by d'^^_^{z), d'^j^{z) ( |6.28|) and 

r(. + i) ' r(2z + i) ^7 20) 

(with Ur taken from ( |6.29D ). The difference equations ( |7.6| ), ( [7.7| ) follow again from 
the standard difference equation for the gamma function. The normalization theorem 
reads in this case. 

Theorem 7.5. Let us assume continuous Hahn parameters with values taken from 
the domain indicated in Section ^.3l such that the recurrence relations of Theorem \6.l\ 
hold. Then one has 



(7.21) 



{pf,pf)A^« _ Af (p'-^ + A) Ai^(p-^ + A) 
(1, 1)a^« Afip^^) A=_^ (p-^) 

(with p'^^ given by (|5.18| ) ). 

For n = 1 the r.h.s. of ( |7.21| ) becomes 

{pf,pf)A^«/{pf,pf)A^« 

with (cf. |[KS|| and recall our normalization in ( |2.18| )) 

{pf,pf)A^«= (7.22) 

(z/o+ + + + 1^1+1- l)i r(z/+ + + Z/+ + z/f + 21) ■ 
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7.4. Jacobi type. For the Jacobi type polynomials the functions Aj-{x) 
of the form in ( [7.5D are determined by d:l,^{z), d:^ _^_{z) ( |6.3(j| ) and 



T{-u + z+l) 
Tiz + l) ■ 



r(-z>o + z + i)r(-t)i + ^ + 1^ 
r(2z + 1) 



(7.23) 

(with z>o, z>i taken from ( |6.37D ). The difference equations (|7!6|) , ( [7.71) follow again 
from the difference equation for the gamma function. The normalization theorem 
reads in this case. 

Theorem 7.6. Let us assume Jacobi parameters with values taken from the domain 
indicated in Section \2.A Then one has 



24|A| 



(1,1)a^ 
(with p"^ given by (|5.21|) }. 
For n = 1 the r.h.s. of ( 7.24| ) becomes 

with (cf. ||AS| , [KS|| and recall our normalization in ( p.21| )) 



(7.24) 



2"'+^a-^/!r(z/o + l/2)r(z/i + l/2) 
{uo + 1/1 + /)/ T{uo + UI + 21 + 1) ' 



(7.25) 
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